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Abstract 

The level set of an elliptic function is a doubly periodic point set in 
C. To obtain a wider spectrum of point sets, we consider, more gen- 
erally, a Riemann surface S immersed in and its sections ("cuts") 
by C. We give S a crystallographic isometry in by defining a fun- 
damental surface element as a conformal map of triangular domains 
and S as its extension by reflections in the triangle edges. Our main 
result concerns the special case of maps of right triangles, with the 
right angle being a regular point of the map. For this class of maps 
we show that only seven Riemann surfaces, when cut, form point sets 
that are discrete in C. Their isometry groups all have a rank- four 
lattice subgroup, but only three of the corresponding point sets are 
doubly periodic in C. The remaining surfaces form quasiperiodic point 
sets closely related to the vertex sets of quasiperiodic tilings. In fact, 
vertex sets of familiar tilings are recovered in all cases by applying 
the construction to a piecewise flat approximation of the correspond- 
ing Riemann surface. The geometry of point sets formed by cuts of 
Riemann surfaces is no less "rigid" than the geometry determined by 
a tiling, and has the distinct advantage in having a regular behavior 
with respect to the complex parameter which specifies the cut. 
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1 Introduction 

Crystallography is concerned with point sets in R" that are discrete and 
distributed more or less uniformly. In "classical" crystallography, periodicity 
was imposed as well, but this restriction is not considered as fundamental 
in the "modern" era. With the discovery of intermetallic quasicrystals |I| 
in the early 1980s, it became clear that there exist aperiodic point sets that 
share a basic property with periodic point sets. This property should really 
be associated with a distribution: in this case, the distribution formed by 
placing a Dirac delta at each point of the set. In these terms, the class of 
aperiodic sets singled out by crystallography is characterized by the property 
that the Fourier transform of the corresponding distribution has support on 
a lattice . The rank of this "Fourier lattice" equals the dimension of space 
for periodic sets, and exceeds it (but is still finite), in the case of quasiperiodic 
sets. The vertex set of the Penrose tiling of the plane is a familiar example 
of a quasiperiodic set. 

The standard construction of quasiperiodic sets A begins by embedding 
M" = F in a larger Euclidean space, M™+" = X x Y. Into X x Y one 
then immerses a smooth m-manifold S that is (z) transversal to Y, and 
(ii) invariant under the action of a lattice A generated by n + m linearly 
independent translations in X x y. Point sets A G Y are obtained as 
sections ("cuts") of S by spaces parallel to Y. More formally, in terms of the 
standard projections 

Tvx-.S^X (1) 

TTY-.S-^Y, (2) 

the section of 5 at s G X is the set 

A{x) = tty o nx^^ix). (3) 

Periodicity or quasiperiodicity of A{x) is determined by the rank of the 
lattice Ay = AdY. Since the generators of A were assumed to be linearly 
independent, rk (Ay) < n. Quasiperiodicity corresponds to rk (Ay) < n, with 
complete absence of periodicity characterized by Ay = {0}. An important 
motivation for constructing quasiperiodic sets, in this context, is the fact 
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that symmetry groups that cannot be reahzed by periodic point sets in M", 
can be reahzed by periodic surfaces in M'""''"'. 

Transversahty and periodicity are relatively mild restrictions on the man- 
ifold S, called the "atomic surface" by physicists. A further restriction, one 
which leads to point sets called "model sets" , is to require that S is the A- 
orbit of a polytope in X. The algorithm which constructs A{x) from such S 
naturally leads to the terminology "window" or "acceptance domain" for the 
corresponding polytopes. Model sets can always be organized into finitely 
many tile shapes, and, because of this simplicity, have dominated the study 
of quasiperiodic sets. 

A different viewpoint on the construction of S, pioneered by Kalugin ^ 
and Katz 0, emphasizes the continuity properties of A{x) with respect to 
X. Consider in more detail the construction of a model set: S = V + A, 
where P C X is a polytope. Now, if x G V + ttxW for some A G A, 
then y = Try (A) G A{x). But now consider what happens when x crosses 
the boundary of V + nx{X)- As x "falls off the edge of the earth", the 
corresponding point y in the point set A{x) disappears. By the same process, 
of course, points can spontaneously appear "out of thin air". To gain control 
over these processes, Kalugin and Katz advocated a restriction on 
V, in relation to A, such that whenever x falls off the edge of one polytope, 
V+7ix{X), it falls within another, say V+itx{X'). This restriction corresponds 
mathematically to the statement that the boundaries of the disconnected 
components of S = V + A can be "glued" together to form a topological 
manifold without boundary. 

In the process of restoring transversahty to the glued complex of polytopes 
one encounters the problems addressed by singularity theory. The map ttx 
should now be a smooth (but not necessarily 1-to-l) map of m-mainfolds. In 
the trivial situation, when ttx has no singularities, S must be diffeomorphic 
to a collection of hyperplanes. This is the situation explored by Levitov |0] 
for point sets in two and three dimensions and various symmetry groups. 

When iS is a generic 2-manifold, we have the classic result of Whitney ^ 
that the stable singularities of smooth maps, such as nx, are folds and cusps, 
having respectively codimension one and two. Because the cusp is always 
accompanied by two folds, the locus of singular values of tcx consists of curves. 
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The space X is thus populated by singular curves such that whenever x 
crosses a curve, a pair of points in A{x) merge and annihilate. One motivation 
for the present work was the desire to eliminate this point-merging singularity 
to the greatest extent possible. 

By giving the 2-manifold S a complex structure, and identifying X with 
the complex plane, we impose additional regularity by insisting that ttx is 
locally holomorphic. The singularities of ttx will then be isolated points. 
A construction that naturally leads to a ttx with this property is to let S 
be (locally) the graph of a holomorphic function, f:X Y. Globally 
this corresponds to a Riemann surface S immersed in C"^ = X x Y and 
having an atlas of compatible charts in X. The other ingredient needed by 
our construction is some way to guarantee that S is invariant with respect 
to a lattice A. We meet this challenge by using conformal maps between 
triangles to define a fundamental graph of S. Schwarz reflections in the 
triangle edges extend this graph and generate the isometry group of S. For 
appropriate choices of triangles, the isometry group has a lattice subgroup 
with the desired properties. 

In the second half of this paper we classify a subset of all Riemann sur- 
faces generated by conformal maps of triangles. This subset is characterized 
by the property that the conformal map is regular at one vertex of the trian- 
gles and that the edges at this vertex make the largest possible angle, 7i/2. 
With the only other restriction being that the corresponding point sets A{x) 
are discrete in Y, one arrives at a set of seven surfaces. Four of these are 
quasiperiodic. The point set obtained from a section of one of them is shown 
in Figure 1. Also shown in Figure 1 is a much studied model set 0: a tiling 
of boats, stars, and jester's-caps (whose vertices coincide with a subset of the 
Penrose-tiling vertex set). The point set determined by the Riemann surface 
can be said to be approximated by the vertex set of the tiling by a systematic 
process that renders the Riemann surface piecewise flat. The other surfaces 
obtained in our partial classification, when flattened, also produce familiar 
tilings (Fig. 3). 
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Figure 1: Point set (small circles) given by the section of a Riemann surface 
compared with a tiling of boats, stars, and jester's caps 0. 
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2 Riemann surfaces generated by conformal 
maps 

2.1 Immersed Riemann surfaces 

We consider Riemann surfaces as analytically continued holomorphic func- 
tions, interpreted geometrically as surfaces immersed in C^. Our treatment 



follows closely the notation and terminology of Ahlfors [10 



Definition 2.1. A function element F = {U, f) consists of a domain U G C 
and a holomorphic function / : U ^ C. 

Definition 2.2. Function elements Fi = {Ui, fi) and F2 = (f/2, /2) are direct 
analytic continuations of each other iS V = Ui H U2 7^ ^ and fi = f2 when 
restricted to V. 

Definition 2.3. The complete, global analytic function determined by func- 
tion element Fq = {Uq, fo) is the maximal collection of function elements 
such that for any Fi E there exists a chain of function elements Fq, . . . ,Fi, 
all in JF, with every link in the chain a direct analytic continuation. 

Up to this point the set of function elements comprising a complete, global 
analytic function JF only possesses the discrete topology, where Fi fl F2 = 
whenever Fi and F2 are distinct elements of JF. By refining this topology we 
can identify J-' with a surface and, ultimately, a Riemann surface. Consider 
a pair of function elements in JF, Fi = (Ui, fi) and F2 = {U2,f2)- In the 
refined topology we define the intersection by 



{.U3, fs) if Fi and F2 are related 
by direct analytic contin- 
uation, (4) 

otherwise, 



where U3 = f/i nf/2 and /s is /i = /2 restricted to U3. It is straightforward to 
check that this defines a valid topology. Moreover, the projection n: T C 
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given by 

n:iU,f)^U, (5) 

provides the complex charts that identify with a Riemann surface. 

Throughout the rest of this paper we will mostly be interested in Riemann 
surfaces immersed in C^. 

Definition 2.4. Let ^ be a complete, global analytic function. The im- 
mersed Riemann surface S corresponding to T is the image of the immersion 
: ^ ^ given by 

([/,/) ^{(a;,/(a;)): a; (6) 

Notation. We denote the first component of by X, the second by Y. 

If we restrict the immersion ^ to a single function element, Fq — {Uq, /o), 
we obtain the graph 

So = {ix,foix))eXxY:xeUo} (7) 

Thus So represents a piece of S and in fact determines all of 5; 5 is connected 
because every pair of function elements in a complete global analytic function 
is related by a chain of direct analytic continuations. S is the completion of 
So- 

Notation. We write [Sq] to denote the completion of the graph Sq. 

Since all subsequent references to "Riemann surface" will be as a surface 
immersed in XxY, wc drop the qualifier "immersed" below. We also omit the 
term "complete" , since the only instances of incomplete surfaces, graphs, will 
always be identified as such. Given a Riemann surface S, we will frequently 
make use of the projections 

TTx-S-^X, (8) 

tty-.S^Y. (9) 

The historical construction of Riemann surfaces we have followed can be 
criticized for its inequivalent treatment of the spaces X and Y. We can 
correct this fault by insisting that the functions / appearing in the function 
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elements [U, /) are not just holomorphic in their respective domains U, but 
conformal (holomorphic with holomorphic inverse). The graph (|^) could 
then be equally written as 



where Vq = fo{Uo). If this "inversion" , or interchange of X with Y, is to work 
for all function elements {U,f), then one must remove all points xq G U, 
where / behaves locally as /(x) — f{xo) = c{x — Xq)*" + ■ ■ • , with m > 1. 
These correspond to branch points of the map vry. Conversely, had we begun 
with the inverted function elements our Riemann surface would have included 
branch points of nx, i-e. points where / is singular. In keeping with tradition 
we augment our definition of a Riemann surface S to include all points (xq, yo) 
where S behaves locally like the algebraic curve (y — yo)"' = c(x — xq)™, where 
m and n are positive integers. 

Definition 2.5. A point (x, G 5 is regular if the corresponding complete, 
global analytic function contains a function element {U, /), with x G f/ and 
/ conformal at x. A point which is not regular is singular. 

2.2 Transformations 

Two transformations of Riemann surfaces will be needed in our discussion 
of symmetry properties. These are defined in terms of their action on the 
spaces X and Y and induce a transformation on Riemann surfaces as subsets 
of X X y . Let (x, y) be a general point in X x Y and define the following 
t r ansf or mat ions : 



Transformation r (for complex constants a, b, c, and d) is the general bilin- 
ear map while a corresponds to Schwarz reflection (componentwise complex 
conjugation). 

Lemma 2.1. If S is a Riemann surface and T is either of the transforma- 
tions T or a, then TS is again a Riemann surface. 



So = {ifo-\y),y)eXxY:yeVo} 



(10) 





(11) 

(12) 
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Proof. Write ^(x, y) — {Txx, Tyy) where Tx and Ty are just maps of the 
complex plane. Since S corresponds to a complete global analytic function J^, 
we need to verify that TS corresponds to some other complete global analytic 
function J-'^. From our definitions we see that J-'t is obtained from JF by sub- 
stituting each function element F = {U, f) E J-' by Ft = {TxU, Ty o f oT^^). 
It is easily checked that Tx is open and Ty o f o T^^ is holomorphic for both 
of the transformations being considered. Thus Ft remains a valid function 
element. One also verifies that the direct analytic continuation relationships 
among function elements are unchanged by these transformations. □ 

Corollary 2.2. If Sq is a graph and T is either of the transformations r or 



Two special transformations are rotations and translations, for which we 
introduce the following notation: 



More generally, transformations T: XxY^XxY which act isomet- 
rically on the spaces X and Y are just the products of Euchdean motions 
in X and Y. Isometries of Euchdean spaces normally include reflections; to 
preserve the structure of the immersed Riemann surface, however, any re- 
flection in X (complex conjugation) must be accompanied by a reflection in 
Y. 

Definition 2.6. The isometries of X xY is the group of transformations 
generated by cr, r(^,0), and t{u,v). 

In what follows we use the term "isometry" only in this sense. Isometries 
which preserve a Riemann surface <S are called isometries of S and form a 
group. The maximal group of isometries is called the isometry group of S. 

Definition 2.7. The group of proper isometries of X x F is the normal 
subgroup of isometries generated by r[6,(l)) and t{u,v). Any element of the 
coset, ag, where g is a. proper isometry, is called a Schwarz reflection. 



a, then [TSq] ^T[Sq]. 



r{e,(l)) = r(e*^ 0; e^"^, 0) 
t{u, v) — r(0, u; 0, v). 



(13) 
(14) 
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2.3 Surfaces generated by conformal maps of triangles 

We now focus on the class of Riemann surfaces determined by graphs which 
solve a purely geometrical problem: the conformal map between two bounded 



triangular regions, P C X and Q CY. The Riemann mapping theorem pH 
asserts there is a three-parameter family of conformal maps f : P Q that 
extend to homeomorphisms of the closures P and Q. To fix these parameters 
we require that the three vertices of P map to the vertices of Q. This defines 
the graph 

P\Q = {{x,f{x)):xeP}, (15) 

and a corresponding Riemann surface [-P|Q]. The closure of P\Q is defined 
analogously and is written P\Q. One of the main benefits of using a conformal 
map of triangles to determine a Riemann surface S is that its isometry group 
can be understood simply in terms of its action on a partition of S into tiles. 

Just as P can be decomposed into an interior P, edges which bound P, 
and vertices which bound each edge, there is a corresponding cell decomposi- 
tion of the graph P\Q. For example, if Pi is one vertex of P, and /(-Pi) = Qi 
is its image in Q, then we use the symbol Pi\Qi to represent the correspond- 
ing vertex of P\Q. Each vertex of P\Q is associated with two angles, a vertex 
angle of P and the corresponding vertex angle of Q. Let the three angle pairs 
be ctj. Pi, i = 1,2, 3. If a, = (3i for all i, then P is similar to Q and / is just 
a linear map. Because the corresponding Riemann surface would be trivial 
(a plane) we exclude this case. It is impossible to have 7^ Pi for just one 
i since then the angle sum could not be vr in both triangles. Thus we must 
have at least two vertices with unequal angles. At these vertices / fails to 
be conformal. Any vertex of P\Q where the corresponding angles in P and 
Q are unequal will be called a singular vertex. The singular vertices of P\Q 
are the only singular points of P\Q. 

The edges of P\Q (associated with each pair of vertices ij = 12, 13,23) 
are effectively the generators of the isometry group of [-PjQ]. By PulQu we 
mean the graph given by the restriction of / to the edge P12 of P with image 
an edge Qu of Q. Consider the triangles P' and Q' obtained from P and Q 
by reflection in these edges. The graph P'\Q', determined by the conformal 
map g: P' —>■ Q', is clearly related to P\Q by an isometry of X xY , a Schwarz 
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reflection which we call ai2- Because au fixes every point of PulQu, we have 
that f{x) = g{x) for all x G Pu- A basic result from complex analysis then 
tells us that the function elements (P, /) and (P', g) are related by analytic 
continuation. Thus [P\Q] = [P'\Q'] = [(ri2{P\Q)] = (ri2[P\Q], by Corollary 
23. 

Definition 2.8. The group G generated by the Schwarz refiections aij which 
fix the three edges of a triangular graph P\Q is called the edge group of P\Q. 
The edge group of P|Q is a subgroup of the isometry group of [P|<5]. 

In order to show that the edge group of a triangular graph is the maximal 
isometry group, we first need to refine the sets on which these groups act. 

Notation. The symbol S corresponds to the Riemann surface S whose sin- 
gular points have been removed. 

Definition 2.9. A real curve of the Riemann surface S is any curve F C 
S, homeomorphic to M, and pointwise invariant with respect to a Schwarz 
refiection. 

Since both ttx'- S ^ X and ny '■ S Y are immersions, the map vry o 
TTx^^ : 7rx(r) —>■ 7ry(r) is an immersion as well. Thus it makes sense to use 
our graph notation, F = 7|5, for real curves, where 7 = vrx(r) and 6 = vry (F). 
A real curve 7|(5 is geometrically no different from the edge of a triangular 
graph; the projections 7 and 6 are always straight lines. Any real curve is 
isometric with the graph of a real analytic function. 

The three real curves which bound the triangular graph P\Q generate a 
topological cell decomposition of P|Q into vertices, edges, and the graph P\Q 
itself. The generators of the edge group, a^j, acting on P\Q, generate three 
closed graphs, each having one edge in common with P\Q. By continuing this 
construction we obtain a cell decomposition of [P\Q] into 2-cells isometric 
with P\Q, 1-cells isometric with one of the edges of P\Q, and points. The 
cell complex as a whole defines a tiling T; the 2-cells by themselves form a 
set of tiles, 7^, and every element of 7^ can be expressed as g{P\Q), where g 
is an element of the edge group, G. 

To show that G is the maximal isometry group we first need to check 
that the tiling T is primitive, that is, there is no refinement of the tiles 7^ by 
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additional real curves within [P\Q] we may have missed. For this it suffices 
to check that there are no real curves within P\Q. Before we can prove this 
statement we need some basic properties of real curves. 

Definition 2.10. A real curve is complete if it is not a proper subset of any 
other real curve. 

Lemma 2.3. The closure in S of a complete real curve 7|5 C S, if bounded, 
has singular endpoints. 

Proof. Without loss of generality let 7 and S lie on the real axes of, respec- 
tively, X and Y. The functions / of the function elements {U,f), which 
represent S locally, will then have power series on the real axis (of X) with 
real coefficients. Since a real power series when analytically continued along 
the real axis continues to be real, we can continue until we encounter 
either a singularity of / or a zero of /' {i.e. a singularity of f~^ on the real 
axis oiY). □ 

The next Lemmas deal with the angles formed by intersecting real curves. 

Definition 2.11. The angle between lines 7 and 7' (in X or Y), denoted 
Z(7, 7'), is the smallest counterclockwise rotation required to make 7 parallel 
to 7'. 

Lemma 2.4. If real curves C S and ^'\5' C S intersect, then Z(7, 7') = 
Z(5,5'). 

Proof. Near the point of intersection S is represented by a function element 
{U, f) where / is conformal. The equality of angles, formed by a pair of lines 
in X and their images by / in y, is simply the geometrical statement that 
/ is conformal. □ 

Lemma 2.5. Only a finite number n > 1 of real curves can intersect at any 
point of a nontrivial Riemann surface and the angle formed by any pair must 
be a multiple o/tt/ti. 
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Proof. Suppose and ^'\5' intersect with angle Z(7, 7') = Z(5, 5') — a> Q 
on a nontrivial Riemann surface S\ for convenience, let (0, 0) be the point of 
intersection. These curves are fixed by Schwarz reflections a and a' respec- 
tively, and a'a = r{2a,2a) is an isometry of S. The neighborhood of the 
point of intersection is the graph 

So^{{x,f{x)):xeU}, (16) 

where U is a neighborhood of the origin in X, and / is conformal at a; = 0. 
The Taylor series for / at the origin has the form 

00 

/(x) = J]a,x^ (17) 

k=l 

where Oi 7^ 0. A short calculation shows 

r{2a, 2a)So = {{x, Ux)) : x G U^}, (18) 

where Ua — e*^"C/ is again a neighborhood of the origin, and 

00 

/e.(x) = J]afee'2"(i-^)a;^ (19) 
fe=i 

Since r{2a,2a) is an isometry, the Taylor series for / and fa must agree, 
term by term. Now if a = ttuj and u is irrational, then uj{l — k) can be an 
integer only for k = 1 (so that e*^"*^^^^) = 1). But this requires = for 
k > 1 which is impossible since S is nontrivial. Thus we may assume lj = p/q 
where p and q are relatively prime positive integers, p < q (since a < tt), and 
di+mq 7^ for some integer m > 0. 

Now let 7" 1 5" be any real curve that intersects 7|(5 at the origin; then 
^(7)7") — ^(<^)<^") — '^{p'/q') by the argument just given, where p' and q' 
are relatively prime positive integers, p' < q'. However, since ai+mq 7^ 0, we 
must have ^-'^'^^ip'/i')'^'} = 1^ or that q' divides the product mq — n. This 
shows that Z{j, 7") is a multiple of tt/ti, for some n > 1. □ 

Clearly any triangular graph with a nontrivial isometry must be "isoceles" 
and fails to be primitive because it can be decomposed into two isometric 
tiles. This is made precise by the following Lemma. 




Figure 2: Diagrams used in the proof of Lemma 
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Lemma 2.6. Let P\Q be a nontrivial triangular graph with trivial isometry 
group, then P\Q contains no real curves. 

Proof. Suppose P\Q contains a real curve and call its completion 'j\5. We 
recall that 7, and its closure in X, 7, are straight lines and 7 cannot have 
an endpoint within P (Lemma p.3| ). The possible geometrical relationships 
between 7 and P are diagrammed in Figure 2. Since P\Q is nontrivial, at 
least two vertices are singular and are shown circled in each diagram. Either 
7 intersects two edges of P, as in cases A and B, or, it intersects an edge 
and the opposite vertex which may be singular (case C) or possibly regular 
(case D). The vertex labels on the diagram refer to our notation for the 
vertex angles and edges. For example, ai and Pi are the angles in P and Q, 
respectively, of vertex 1; P^lQu is the edge (real curve) bounded by vertices 
1 and 2, etc. 



Case A is easily disposed of using Lemma 



Z(7, P12) = Z(7, P13) + «i = Z(5, Q13) + «i = Q12) -Pi + ai (20) 

= Z(7,Pi2)-/?i + «i. (21) 

This is impossible because vertex 1 is singular (ai 7^ Pi). 

By using Schwarz reflection to imply the existence of additional real 
curves, the remaining cases either reduce to case A or imply the existence of 
a singularity within P\Q ot one of its edges — neither of which is possible. 

First consider case B. Let 7 intersect P13 at xi and P23 at X2, forming 
angles 61 and 62 (see Fig. 2). Any other complete real curve with projection 
7' which intersects Xi makes a finite angle with 7 by Lemma |2.5| . Thus 
we may assume the angles 61 and 62 are the smallest possible (for a 7 that 
intersects both P13 and P23). Since one of 9i and 62 must be greater than 
n/2, we assume without loss of generality it is 9i. If we now reflect 7 in P13 
we obtain a real curve with projection 7' such that 7' intersects P13 but not 
P23. Thus case B always reduces to cases A or C. 

In case C we consider the sequence of real curves with projections 7^, 
where 70 = P23, 71 = 7, and 7^+1 is the image of 7^-1 under reflection in 
7fc. Let 6k be the angle formed at vertex 2 in P by 7^. Clearly for some 
k we arrive at a 7' = 7^ such that 9 = 9k < 0:2/2 (see Fig. 2). This 
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leads to three subcases: Ci, where Z(7', P13) < 7r/2, C2, where 7' and P13 
are perpendicular, and C3, where Z(7',Pi3) > n/2. In case C2, 9 < a^jl 
since otherwise P\Q would have a nontrivial isometry (reflection in 7'|5')- 
All three subcases immediately lead to contradictions. In Ci, reflecting 7' 
in P13 presents us with a 7" satisfying case A. In C3, the image of vertex 1 
under reflection in 7' implies a singularity within P; in C2 the same reflection 
implies a singularity on P13. 

Case D: we either have Z(7,Pi2) = 7r/2, case Di, or Z(7,Pi2) 7^ 7r/2, 
case Since P has no nontrivial isometry, a reflection in 7 in case D\ 
would place the image of either vertex 1 or 2 (both singular) somewhere on 
Pi2. In D21 a Schwarz reflection of 7 leads to case A. □ 

Theorem 2.7. Let P\Q he a nontrivial triangular graph with trivial isometry 
group, then the maximal isometry group of [P\Q] is the edge group of P\Q. 

Proof. We use the real curves to decompose [P\Q] into a set of tiles (2-cells) 
T2. Lemma tells us that P\Q G T2, so that T2 = G{P\Q), where G is the 
edge group of P\Q. On the other hand, if h is an isometry of [P|Q], then 
h{P\Q) = P'\Q' e r2, where P'\Q' = g{P\Q) for some g e G. But since 
P\Q has no nontrivial isometries, the map h~^g: P\Q P\Q must be the 
identity and h = g. □ 

We conclude this section with a formula for the topological genus of a Rie- 
mann surface [P\Q] compactifled by the translation subgroup of its isometry 
group, the lattice group A of [P|<5]. 

Definition 2.12. The vertex groups Gi, {i = 1,2,3), of a triangular graph 
P\Q, are the subgroups of the edge group of P|Q generated by the adjacent 
edges of, respectively, the three vertices of P\Q. 

Theorem 2.8. Let P\Q be a nontrivial triangular graph with trivial isometry 
group. Let G be the isometry group of [P\Q], A its lattice group, and Gi, 
{i = 1, 2, 3), the three vertex groups of P\Q. If \G/A\ is finite, the genus g of 
the surface [P\Q]/A satisfies 




(22) 
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Proof. If I G/ A I is finite we can view [P|Q]/A as a finite cell complex. We can 
relate the number of 0-cells, Nq, and the number of 1-cells, A^^i, in this complex 
to the number of 2-cells, N2. Since every 2-cell is bounded by three 1-cells, 
each of which bounds exactly one other 2-cell, A^i = {3/2)N2. Similarly, the 
boundary of each 2-cell contains three 0-cells (the vertices i = 1,2,3), each 
of which belongs to the boundary of a number of 2-cells equal to the order of 
the corresponding vertex group, \Gi\. Thus A^o = (X^Li \Gi\~^)N2. Finally, 
since P\Q has no nontrivial isometry, and G/A acts transitively on the 2- 
cells of [P\Q]/A, N2 = \G/A\. The result follows from Euler's formula, 
2 ~2g = N2- Ni + Nq. □ 



2.4 Discreteness and uniformity 

The whole point of immersing a Riemann surface S in X xY = C"^ is that by 
forming sections of S, i.e. intersections with X = C, one obtains patterns of 
points. A very primitive property of a point set, normally taken for granted 
in crystallography, is discreteness. 

Definition 2.13. The section of the Riemann surface 5 at x is the set 

Aix) =7lYOTlx~\x). (23) 

The basic property of a holomorphic function, that its zeros form a dis- 
crete set, translates to the statement that the preimages nx~^{x) are discrete 
in a Riemann surface S. We use the stronger property that ttx~^{x) is dis- 
crete in X X y to define a discrete Riemann surface. This is equivalent to 
the following statement about sections: 

Definition 2.14. A Riemann surface S is discrete if its sections A{x) are 
discrete in Y for every x G X. 

All the statements we can make about discreteness of a Riemann surface 
S hinge upon properties of the lattice group of S, A. One property is the 
rank, rk(A), given by the cardinality of the generators of A. The orbit of 
the origin of X x y, A(0,0), is called a lattice and is also represented by 
the symbol A. When rk(A) = 4, a second property is the determinant of the 
lattice, det A. If det A > 0, the four generators of A are linearly independent 
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(as vectors m X xY); if det A = the generators are linearly dependent and 
A (as a lattice) is not discrete in X xY. A lattice with rk(A) > 4 is never 
discrete in X x Y. 

Notation. The standard measure for a set A is written \A\. If A is a region 
in C then l^l is its area; if A is a set of points, then 1^41 is its cardinality. 
Finally, if A is a rank 4 lattice, then |A| = Vdet A is the volume in X x F of 
its fundamental region. 

The following Lemma provides a necessary condition for discreteness: 

Lemma 2.9. The lattice group A, of a discrete, nontrivial Riemann surface 
S, is discrete (as a lattice) in X xY and in particular, rkA < 4. 

Proof. If A is not discrete we can find a sequence of t{u,v) G A such that 
both —>■ and \\v\\ —>■ 0. Let (a;o,2/o) be a regular point of <S, then 
Ho e A{xo). Near {xo,yo) we can represent S by the graph 



where C/ is a neighborhood of Xq in X, f is conformal in U, and /(xq) = yo- 
Since t{u,v) is an isometry, t{u,v){U\V) C S. Prom 



we see that f{xQ — u)-\-v e ^(xq), since Xq — u G U as \\u\\ can be arbitrarily 
small. Since S is discrete, yo is isolated in Y and there must be a subsequence 
t{u',v') such that f{xo — u') + v' = f{xo). If, within the sequence t{u',v'), 
there is a subsequence t{u", v") with u" = 0, then v" = f{xo)—f{xo—u") = 
and we have a contradiction. Thus there must be a subsequence with u" ^ 0. 
Since / is conformal at Xq, 



U\V^{ix,fix)):xeU}, 



(24) 



t{u,v){U\V) = {{x + u,f{x)+v): X e U}, 



(25) 



lim 



/(^o)-/(^o-^^") 



/'(^o) 



(26) 



u' 







U' 



.11 



(«",i;")^(0,0) V! 



II 



(27) 



But the second limit, above, is independent of Xq so we are forced to conclude 
that /' is constant. This is impossible because S is nontrivial. □ 
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The point sets studied in crystallography normally are Delone sets and 
have the property of being uniformly discrete ^ . For a point set in M" this 
means there exists a real number r > such that a spherical neighborhood of 
radius r about any point of the set contains no other point of the set. For the 
sets A{x) generated by Riemann surfaces this property is clearly too strong: 
it is violated whenever x is near a branch point of ttx- We therefore adopt a 
weaker form of this property which is nevertheless stronger than discreteness 
and useful in establishing the existence of the density. 

Definition 2.15. Let Br{y) C F be an open disk of radius r centered at 
y. A Riemann surface S, with sections A{x), is finitely discrete if for some 
r > 0, |^(a;) fl Br{y)\ is uniformly bounded above for all x G X and y 

Since a disk of radius r' can always be covered by finitely many disks 
of radius r, the finitely discrete property holds for any r' once it has been 
established for a particular r. 

We now introduce the class of Riemann surfaces which is the focus of this 
study. 

Definition 2.16. A Riemann surface is crystallographic if its lattice group 
A has rank 4 and |A| > 0. 

Definition 2.17. Any isometry g of X x Y can uniquely be expressed in 
the form g = goX, where go fixes the origin and A is a translation. The 
derived point group of the isometry group G, ip{G), is the image of G by the 
homomorphism if): g y-^ g^. Since Ker^ = A, the lattice group of G, we have 
the isomorphism ip{G) ~ G/A. 

It is important to remember that ipiG) need not be a subgroup of G\ nev- 
ertheless, the lattice group of G is always left invariant by ip{G), just as it 
is invariant within G. Furthermore, if a Riemann surface is crystallographic, 
then the action (by conjugation) of '>p{G) on its lattice group A is a faithful 
representation of G/A. Since the isometry group of a (finite rank) lattice is 
finite, we have that G/A, for a crystallographic Riemann surface, is always 
finite. 

Lemma 2.10. A Riemann surface determined by a triangular graph P\Q, if 
crystallographic, is finitely discrete. 
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Proof. Let Br{x,y) C X x y be an open ball of radius r centered at an 
arbitrary point {x,y). Consider the piece of the Riemann surface within 
this ball, Sb = [P\Q] H Br{x,y). and the projection Tlx' Sb — ^ X. [P\Q] is 
finitely discrete if there is a uniform upper bound on the number of preimages 

[P\Q] is covered by the orbit of closed graphs, G{P\Q), where G is the 
edge group of P\Q. Let A be the lattice group of G, then G{P\Q) is the union 
of cosets, Hi{P\Q), i — 1, . . . ,N , where N — \G/h\ is finite because [P\Q\ 
is crystallographic. Again, because [P\Q\ is crystallographic, all but finitely 
many graphs in Hi{P\Q) have empty intersection with a ball of radius r, in 
particular, Br{x,y). Thus we have a bound (independent of x and y) on the 
number of graphs in G{P\Q) which intersect Br{x,y). But a graph can have 
at most one preimage of Tlx', hence |7r^^(a;)| is uniformly bounded above. □ 

The sections A{x) of a crystallographic Riemann surface also possess a 
uniformity with respect to the parameter x. Our handle on this property is 
provided, in part, by the smooth behavior of A{x) with x. Before we can 
proceed, however, we need to be aware of two point sets in X which create 
problems: branch points (of nx) and crossing points. 

Definition 2.18. A point {x,y) G 5 is a self-intersection point if in the 
description of 5 as a complete global analytic function there exist function 
elements {U, f) and {V,g), such that xeUHV, fy^gmUnV, and 
f{x) = g{x) = y. The point x & X is called a crossing point. 

Lemma 2.11. A Riemann surface determined by a triangular graph P\Q has 
countably many self-intersection points and nx has countably many branch 
points. 

Proof. [P\Q] is the union of countably many closed graphs Pi\Qi given by the 
orbit of P|Q under the action of the edge group. Since each graph has at most 
three singular points, ttx has countably many branch points. If there were 
uncountably many self-intersection points then uncountably many must arise 
from one pair of distinct graphs, say Pi\Qi and Pj\Qj. Let fi and fj be the 
corresponding conformal maps; then fi{x) — fj{x) would have uncountably 
many solutions x G Pi H Pj. Thus either fi — fj, a contradiction, or the 
zeroes of fi — fj would not be isolated, another impossibility. □ 
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Lemma 2.12. Let [P\Q] be crystallographic, Xc (Z X its crossing points, 
and Xf, C X the branch points of Tlx i then for any pair x,x' & X \ (XhUXc), 
there exists a bisection of sections of[P\Q], 

A{x)^ A{x'), (28) 

such that \\y — '^{y)\\ is uniformly bounded above for y e A{x). 

Proof. We arrive at ^ by composing bijections 

•^V. A{x)^A{x"), (29) 
•^2: A{x")^ A{x'), (30) 

such that ||y — ^i(i/)|| and \\y" — ^2(1/") II (correspondingly) uniformly 
bounded. The Lemma then follows by application of the triangle inequality. 

Since [P\Q] is crystallographic, we can partition X xY into translates 
of a bounded fundamental region, V^(0), of its lattice A. Thus for any pair 
x,x' & X we can write 

(x,0) e V(0) + A, (31) 
(x', 0) e V(0) + A', (32) 

where A, A' e A. Consider the point 

{x\ y") = (x, 0) + A' - A e V{^) + A'. (33) 

Since 

{x" - x\ y") = (x", y") - (x', 0) e ^(0) - ^(0), (34) 
both ||a;" — x'll and \\y"\\ have upper bounds independent of x and x'. Now 

A{x)^{yeY:{x,y)e[P\Q]}, (35) 

and 

A{x) + y" = {y'eY: (x, y' - y") e [P\Q]}. (36) 

But (x,y'-y") e [P\Q] iS (x,y' -y") + X" e [P\Q], where A" e A. Choosing 

X" = X' -X = {x" -x,y"), (37) 
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we obtain 

A{x)+y" = {y' e Y: {x",y') G [P\Q]} = A{x"). (38) 

As our first bijection we take the translation = y + y", where \\y"\\ 

is uniformly bounded from above. The point of this intermediate step is 
that for \E'2 we need consider only pairs of sections with bounded separation 

In constructing \E'2 we avoid branch points and crossing points. Since 
X G X \ (Xb U XJ, equation implies x" eX\ {X^ U XJ. Let 7 : [0, 1] 
X \ (Xfe U Xc) be a smooth rectifiable curve with 7(0) = x" and 7(1) = x'. 
To show that 7 exists we recall that Xf, U X^ is countable. We can then find 
7 in the uncountable family of circular arcs with endpoints x" and x', since 
each point of Xb U Xc can eliminate at most one arc. 

The curve 7(t) generates a homotopy of the sections A{x") and A{x'). 
At each point y" G A{x"), 7(t) is lifted to a unique curve 7(t)|5(t) C [P\Q] 
with endpoint (7(0), 5(0)) = {x",y") and we define our second bijection by 
^2(1/") = (^(l)- To finish the proof we need to show that \\S{1) — S{0)\\ is 
uniformly bounded. 

Let P be the closed subset of P that is a suitably small distance r or 
greater from any of its vertices that are branch points of ttx- Let P\Q C P\Q 
be the corresponding graph. The orbit under the edge group, S = G{P\Q) C 
[P|Q], is a Riemann surface from which all the points of ramification (of the 
map Tlx) have been "cut out". The complement, <S = [P\Q]\S, is the disjoint 
union of the branched neighborhoods of all the points of ramification. It is 
possible to find curves 7(t), such as the circular arcs considered above, where 
the branched neighborhoods Si visited by 7(t)|(5(t) are visited only once, for 
t E Ti d [0,1]. Also, because we can bound the length L of 7 and there 
is a minimum distance between branch points (on the branched covering of 
X), the number of such subintervals Tj is bounded, i.e. i = 1, . . . , N. The 
bounds on N and L are uniform bounds, independent of the points x and x'. 

Two additional bounds are needed before we can proceed to bound ||5(1) — 
5(0) II . The first is an upper bound D on the diameter of the projection 
of a branched neighborhood, 7ry(iSi). This follows from the fact that Si 
is isometric with the branched neighborhood Sj of a vertex of P\Q, and 
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Sj C Gj{P\Q), where Gj is the corresponding vertex group. Clearly the 
maximum diameter of 7ry(Gj(P|(5)) is bounded because Q is bounded. 

The map f : P —>■ Q (which defines P\Q), when restricted to P is confor- 
mal and ||/'|| has a maximum value, /i, since P is closed. This means that if 
-f{t)\6{t) G P|g, then 



dS 



dt 



f'il) 



d'y 
H 



dj 



dt 



(39) 



Because S is generated from P\Q hj the action of G, this bounds applies 
globally, for 7(t)|5(t) G S. 

We are now ready to complete the proof: 



r d6 


< 


r d6 




r d6 


J[o,i] dt 


7u^,T, dt 


+ 


-'[0,1]\U^,T, dt 



11^(1) -5(0)11 

For each piece of the curve in a branched neighborhood we have 

< D, 



dS 
'dt 



while in the complement (S), 



d6 
'dt 



< 



d6 



dt 



1 


(^7 


/[0,l]\u£iT, 


dt 



Inequality (^OD thus becomes 

115(1) -5 



<ND + fiL. 



(40) 



(41) 



(42) 

(43) 

(44) 

(45) 
□ 



For discrete Riemann surfaces with sufficiently uniform sections A{x), 
one can define their density. 
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Definition 2.19. Let -Bij(O) G Y he a. disk of radius R centered at tlie 
origin. The limit 

if it exists and is finite, is tlie density of A{x). 

With the aid of Lemmas |2.1CI| and p.l2| we can show, that for a crystallo- 
graphic Riemann surface generated by a triangular graph, p{x) exists and is 
(essentially) independent of x. 

Notation. The standard volume form in X is ux = dx A dx, its puUback on 
a Riemann surface S is written TTx*i^x- 

Theorem 2.13. If[P\Q] is crystallographic with lattice group A, its sections 
A{x) have density 

P = Tn ^x*uJx, (47) 

J[P\Q]/A 

independent of x, provided x is not a crossing point or a branch point of tix ■ 
If G is the edge group of P\Q, then 

P^». (48) 

Proof. Notation. The expression c = O (l/R) indicates there exist constants 
Ci and C2 (independent of R) such that for sufficiently large R, Ci/R < c < 

C2/R. 

Let -B_r(0) C y be a disk of radius R centered at the origin and let 

Nr{x) = \Bn{0)nA{x)\. (49) 

We first obtain a bound on the difference, Nr{x) — Nji{x'), when neither x 
nor x' is a crossing point or a branch point of ttx- By Lemma |^.12| there 
exists a bijection A{x) A{x') such that if ?/ G -B_r(0) fl A{x), then 
\E'(?/) G i?/i;+rf(0) n A{x'), where d > is a constant independent of R, x, and 
x'. This shows 

NRix) <\Bn+diO)nAix')\ (50) 
= Nn{x) + \{Bn+d{0) \ 5/?(0)) n A{x')\. (51) 
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We can cover the annulus -Bij,+d(0) \ -Bh(O) by Mr disks Br{y') of a fixed 
radius r > 0, where, for sufficiently large R, Mr < niR and m is a constant 
independent of R. By Lemma |2.f 0| , \B,f{y') n^(a;')| < n, where n is indepen- 
dent of x' and y'. Thus Nr{x) — Nr{x') < mnR. Combining this bound with 
the bound obtained by interchanging x and x', we arrive at the statement 



We now introduce a disk Cr^O) C X and consider the region Wr = 
Cr{0) X Br{0) C X X Y. Since the set of branch points and crossing points 
is countable and has zero measure in X, and ttx is otherwise smooth, 



Because [P\Q] is crystallographic, we can partition X x Y into translates of 
a bounded fundamental region of its lattice, ^(0). Let tta : [P\Q] [P\Q]/^ 
be the standard projection on the quotient. On [P\Q] H ^(0) the map tta is 
1-to-l and 



This defines p, which we can make positive by appropriate choice of orienta- 
tion on [P\Q]. 

Turning now to the region Wr, there is a maximal subset A_ C A such 
that A_ + 1^(0) C Wr and a smallest subset A+ C A such that Wr C 
A+ + V{0). If A G A+ \ A_, then Vx = iX + V{0)) HWr is a proper subset of 
a fundamental region and 



NRix) - NRix') = \BR{0)\Oil/R). 



(52) 




(53) 




(54) 




(55) 



From this it follows that 




(56) 



and, from straightforward estimates of A_|_ and A_, we conclude 




(57) 
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The projection 7rx([-P|<5] riWi^) covers the disk Cr{0) multiple times, the 
multiplicity at the point x G C/j(0) being the number Nji{x) defined above. 
Thus 



7Tx([P\Q]nWR) 



Nr{x)ujx- 



(58) 



Cfl(O) 



We can again neglect the countable set of branch points X}, and crossing 
points Xc to argue, for xq E X \ (X^ U X^) fixed. 



Nr(x)ujx 



Cr(0) 



Nr{xo)uJx+ / {Nr{x) - NR{xo))iUx 
Cr{0) JCr(0) 

NR{x,)\CRm + |CR(0)||5«(0)|O(l/i?), (59) 



where in the last step we used (0). Combining (p3D, 
using \Wr\ = \Br{0)\\Cr{0)\, we obtain 



, ([59|), and 



Nr{xo) 

\BRm 



p(l + 0(l/i?)). 



and thus 



Nr{xo) 
hm — 7 = p. 

R^oo \B 



(60) 



(61) 



To evaluate p from (|5^) , we regard [P\Q]/A as |G/A| equivalence classes of 
tiles, all isometric to P\Q. The result (^8]) follows because the integral of the 
form 7rx*uJx over P\Q is just the volume of ttx{P\Q) = P in X. □ 

Formula (|^) for the density was introduced by Kalugin p to extend the 
notion of stoichiometry to quasicrystals. Because the 2-form 'n'x*(^x is closed, 
this formula gives the same density for 2-manifolds homologous in the torus 
{X xY)/A. One must remember, however, that this homology invariant only 
corresponds to the true density when the map irx is orientation preserving 
(see equation (^)). 
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3 Classification of discrete Riemann surfaces 
generated by conformal maps of right tri- 
angles 

3.1 Conformal maps of right triangles 

The simplest nontrivial conformal maps of triangles, f: P ^ Q, are those 
where one of the vertices of the corresponding graph P\Q is regular. The 



edges of P\Q adjacent to this vertex are real curves, and by Lemma |2.5| , 
belong to a set of n > 1 real curves intersecting at the same vertex with 
minimum angle n/n. This suggests that among those maps with one regular 
vertex, the simplest case is n = 2, i.e. the conformal map of right triangles 
with the right angle being the regular vertex. 

Without loss of generality, we give P and Q a standard scale, position, 
and angular orientation as specified by the vertices Xi\yi of the corresponding 
graph P\Q: 

xi\yi = 0|0 

2^21^2 = cos a I cos /3 (62) 

where < a < n /2 and < /3 < 7r/2 are two free real parameters. Below we 
frequently use the abbreviations a = cos a, b = cos (3. P and Q have angles 
a and /3, respectively, at vertex 1, and the corresponding complementary 
angles at vertex 3. A nontrivial graph P\Q has a ^ (3 with vertices 1 and 3 
singular; vertex 2 is regular. 

The Schwarz-Christoffel formula |1^ gives the conformal map / as the 



composition f = ho g ^ where g and h map the upper half plane of Z 
conformally onto, respectively, P and Q. Explicitly: 



g{z) =A z--\l - z)-2dz, (63) 
Jo 

h{z)=B I z^-\l~z)-hz. (64) 







In both ( |55D and (0) the branches of the fractional powers are chosen so that 
the integrands are real and positive for z G (0, 1). The normalization factors 
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A and B are positive real numbers determined by the conditions g{l) = a 
and h{l) = b. Further properties of g and h are easily checked, in particular, 

g{oo) = 6*°" and h{oo) = e^^ . 



3.2 Isometry groups 

Notation. Let K he a. set of elements of a group G, and let A; G G be some 
element. We denote by {K) the subgroup generated by the elements of 
and by {k}G the conjugacy class of k in G. 

In the case of right triangles, P\Q can have a nontrivial isometry only if 
a = 7r/4 and (3 = n/A. But this makes P\Q trivial. From Theorem we 
know that for nontrivial P\Q the isometry group G is just the edge group 
generated by the three Schwarz reflections: 

0-12 = cr, (65) 

ai3 = r(2a,2/5)(T, (66) 

(T23 = t(2a, 26)r(7r, 7r)(T. (67) 

From the isometry group 

G = {ai2, ai3, (T23) = (a, r(2a, 2/3), t(2a, 26)r(7r, tt)), (68) 

we wish to extract the lattice group A. Helpful in this enterprise are the 
vertex group 

Gi = {a,r{2a,2P)), (69) 

and its cychc subgroup, 

R= {r{2a,2(3)). (70) 

Sets of translations invariant with respect to Gi, or stars, play a central role 
in the construction of A. In what follows we will need two stars: 

S = {t(2a,26)}Gi, (71) 
= r(7r, tt) S r(7r, vr) = {t(-2a, -26)}g,. (72) 



The main result is contained in the following Lemma: 
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Lemma 3.1. Let G he the isometry group of the Riemann surface [P\Q] 
generated by the conformal map of the right triangles specified in (|62|) and let 
R be defined by Gi by (H), S andJ:-^ by and (^. Ifr{7i,n) e R, 
then G has lattice group A = (S) and G = AGi; otherwise, G has lattice 
group A = (SS'i) and G = AGi U t(2a, 2b)r{n, n)AGi. 

Proof. First consider the case r(7r, vr) G R; then 

G = {a,r{2a,2p),t{2a,2b)). (73) 

Now consider the group H = AGi, where A = (S) is normal in H. Clearly 
H G G. Moreover, one easily verifies gH = Hg = H, where g is any of the 
three generators of G. These two facts together show G = if; A is clearly 
the lattice group of G. 

Next consider the case r(7r, vr) ^ R. For the generators of G we must 
now use Consider the group G = AGi, where A = (SS^^) is normal 

in G. Clearly G G G. In contrast to the previous case, we can now only 
verify that gGg^^ = G, where g is any of the three generators in (|68D. Thus 
G is normal in G. G has index at most two, since multiplication of G by 
the generators of G produces at most two, possibly distinct, cosets: G and 
G' = t{2a, 2b)r{TT, 7r)G. But if G' = G, then we would have some A G A and 
some gi G Gi such that \gi = t{2a,2b)r{7i,7!-), or gi = X~H{2a,2b)r{7i,ii). 
Since gi fixes the origin, X~^t{2a, 2b) must be the trivial translation and 
gi = r(7r,7r). This contradicts our assumption r(7r,'7r) ^ R and we conclude 
that G and G' are distinct. Let A be the lattice group of G. Clearly A C A. 
Now suppose A G A but A ^ A. Since then A ^ G, we must have A G G', 
that is, A = t{2a, 26)r(7r, n)gi\ for some gi G Gi and A G A. But this 
implies r(7r,7r)(7i = t(— 2a, — 26)AA^^, a translation, and we arrive at the 
contradiction r{'K^'K)gi = 1. Thus A = A. □ 

3.3 The discreteness restriction 

The requirement that [P\Q] is discrete places strong constraints on the angles 
a and (3 of the right triangles P and Q. 

Lemma 3.2. If [P\Q] is discrete, then a and j3 are rational multiples of it. 
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Proof. By Lemma 2^, [P\Q] has a discrete lattice, i.e. the orbit A ■ (0,0) is 



discrete in XxY. For right triangles. Lemma ^]T| gives us A = (S) if r (tt, tt) G 
R, A = otherwise. Thus discreteness of [P\Q] implies discreteness of 

the star S ■ (0, 0) in X x Y. Since S = {t(2a, 26)}gi = {^(2a, 26)}ij, S ■ (0, 0) 
is just the orbit of t{2a, 2b) ■ (0, 0) under action of the group R generated by 
r(2a,2/5). Clearly S ■ (0,0) lies in a 2-torus Si x Si embedded in X x F. 
Since S ■ (0, 0) is discrete, there is a disjoint union of neighborhoods, each 
containing just one element of S ■ (0, 0). Moreover, since R is an isometry of 
XxY and acts transitively on S ■ (0, 0), there is a uniform lower bound on 
the volumes of these neighborhoods. This implies the existence of disjoint 
neighborhoods in SixSi, again with a uniform lower bound on their measure. 
Since 5*1 x Si has finite measure, this is only possible if R has finite order. □ 

Notation. Given positive integers m and n, GCD (m, n) is their greatest com- 
mon divisor, LCM (m, n) their least common multiple. 

Lemma ^]2| allows us to write a = ir^i/k), j3 = n{j/l), where i, j, k 
and / are positive integers and GCD {i,k) = GCD (j, /) = 1. Since R = 
{r{2TT{i/k),2n{j/l))), we identify n = LCM{k,l) as the order of R. A more 
convenient parameterization is given by 

a = TT — , P = TT — , (74) 

n n 

where p = i{n/k), q = j{n/l) are positive integers with no common factors 
that are also factors of n. 
The transformation 

t(sin a, sin [3) r(— 7r/2, — vr/2) t(— cosa, — cos 13) a, (75) 

has the effect of replacing the angles a and /3bya' = 7r/2— a and (3' = 7r/2— /3 
in the definition of the triangles P and Q. Without loss of generality we may 
therefore take the smallest of the four angles a, a' , (3 and (3' , and rename 
this angle a. Transformation (|75D, as well as an interchange of the spaces X 
and y, will then give all other cases of triangles P and Q (in their standard 
position, orientation and scale). From < a < a' we obtain 

TL 

0<p<-. (76) 
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The inequalities a < (3 and a < (3\ plus the condition to avoid triviality, 
a ^ [3, then give 



n 

p<q< ^-P- 

Together, inequalities ( [75| ) and (|7^) have no solution unless 

n > 6. 

Since the vertex group Gi is generated by R and a, 

\Gi\ = 2n. 



(77) 



(78) 



(79) 

We can use transformation (ffSD, which has the effect of interchanging the 
angles at vertices 1 and 3, to compute the order of G3. Let n', p' and q' be 
the integers parameterizing the angles at vertex 3 (as in ([7^) ); then 

2n 



n 



GCD{2n,n-2p,n-2q) 



(80) 



n n' 2 

i + ^ = 1 
n n' 2' 



and. 



\G, 



2n. 



(81) 



Finally, since the group of the regular vertex is generated by two real curves 
which intersect at right angles. 



I Go 



4. 



From (3S) we find 

7/.(G) = (Gi,r(7r,7r)) = (a,i?,r(7r,7r)) 



(82) 



(83) 



for the derived point group of G. Recognizing n = \R\ as the order of an 
element of the isometry group of the lattice A, we can use the following 
theorem of Senechal |TI|] and Hiller []T5|, and the requirement of discreteness, 
to bound n from above. 
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Theorem 3.3 (Senechal [0, Hiller [0). Let N{n) be the smallest in- 
teger such that the group GL(A^(n),Z) has an element of order n, then 

N{n) = J2 ^(P^^ (84) 

where pi^^ ■ ■ ■ is the prime factorization of n and (f){k) is Euler's totient 
function: the number of positive integers less than and relatively prime to k 
(note: the prime 2 is included in the sum only if its exponent is greater than 

V- 



Lemma 3.4. // a Riemann surface [P\Q] generated by the conformal map 
of right triangles P and Q is discrete, then the angles of P and Q are given 
by (fflD and either n < Q, or n = 8,10 or 12. 



Proof. Let G be the isometry group of [-P|Q], 'ipiG) its derived point group, 
and A its lattice group. Since conjugation hj g & i'iG) leaves A invariant, 
consider the automorphisms $g : A — A given by $g(A) = gXg^^. If A has 
rank N, then the homomorphism \E': ip{G) Aut(A), where "^{g) = 
induces a representation of iplG) by integral N x N matrices of determinant 
±L In fact, is an isomorphism since for any of the generators g of ip{G) 
we can easily find a A G A which is not fixed by $g (it suffices to look 
within the star S or, if r(7r, vr) ^ R, SE~^). Since R C ipiG) has order 



n, there must be an element of order n in GL (A^, Z). By Theorem 3^ we 
must have N > N{n). On the other hand, if [P\Q] is discrete, then A must 
be discrete (as a lattice in X x y) which is possible only if < 4. Since 
0(p"*) = p™^^(]9 — 1), we need only consider the values of for powers of 
small primes: 0(4) = 2, 0(8) = 4, 0(3) = 2, 0(5) = 4 (all other primes and 
higher powers yield values greater than 4). From these facts we obtain just 
the values of n given in the statement of the Lemma. □ 



With Lemma [3^ and inequalities (|76D and (0), the set of possible com- 



binations of n, p and q is already finite. Several of these combinations can be 
eliminated by the following Lemma which provides a lower bound on rk(A) 
when either p or g is a nontrivial divisor of n. 
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Lemma 3.5. Let A be the lattice group of a discrete Riemann surface [P\Q] 
with P and Q defined by parameters n, p and q ([74|),- then if p > 1 and p 
divides n, rk(A) > (f){n/p) + (t){n/ GCD (n, q)) (and the same statement with 
p and q interchanged). 

Proof. Suppose p divides n and n/p = d > 1. We may assume that p 
does not divide q since otherwise n, p and q would have p as a common 
divisor. Let r = r(2a,2/5) = r{27i/d,2n{q/n)); then R = (r) and Rp = {r'^) 
is a subgroup of R of order p. By looking in S (or SE~^), we can find a 
translation to = t{u,v) G A such that m 7^ and v ^ 0. Now consider the 
two products of translations 

t^ = t{uuv^) = Y[{stos~^), (85) 

d 

t2 = t{u2,V2) = YliSktoSk~^), (86) 
k=l 

where Sk is any element of the coset r'^Rp, and ^2 depends on the particular 
choice of coset elements. Evaluating the products we find 

Ml = pu, (87) 

\m=l 



Equation ( pH] ) implies e '^'^^^i/p^vi = Vi, and, since p does not divide q, we 
conclude = 0. Thus ti = t{pu,0). Similarly, we find 



U2 



J2 e'"'^'^/'^) ]u = 0. (89) 



On the other hand, V2 is changed just by making a different choice for one 
coset element (again because p does not divide q). Thus we can always 
make a choice such that ^2 = ^(0, V2), where V2 7^ 0. 

Now {{ti}ji) = Ax is a lattice in X isomorphic to the cyclotomic lattice 
Z[e2^*/'^], while ({^2}/?) = Ay is a lattice in Y isomorphic to Zfe^^*^^/")]. Since 
A D Ax Ay, rk(A) > rk(Ax) +rk(Ay). The statement of the Lemma follows 
from the well known formula for the rank of a cyclotomic lattice. □ 
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As an example of the application of Lemma p.5| , consider the case n = 8, 
p = 1 and q = 2. For these numbers Lemma p75| gives rk(A) > 0(4)+0(8) = 6, 
and the corresponding Riemann surface would not be discrete. Together, 
Lemma |3.4| , Lemma |3.5| and the inequalities ( [f6D and ( ffTl) limit the set of 
possibilities for n, p, and q to the combinations {n,p,q) = (6, 1,2), (8, 1,3), 
(10,1,3), (10,1,4), (12,1,5), and (12,2,3). All other combinations (con- 
sistent with the inequalities) correspond to Riemann surfaces whose lattice 
groups have ranks exceeding 4 and therefore cannot be discrete. We have 
no general method to settle the discreteness of these remaining candidates 
and therefore considered them case by case. The procedure was to obtain 
generators Cj of the finite sets S, or if r(7r,7r) ^ R, since these gener- 

ate the lattice group A. The computations were performed in Mathematica 
using the function LatticeReduce, which implements the Lenstra-Lenstra- 
Lovasz (LLL) lattice reduction algorithm. With the exception of the case 
{n,p,q) = (10,1,4), which was found to have rank 8, all others proved to 
have rank 4 and nonvanishing determinant. To these 5 cases of crystallo- 
graphic Riemann surfaces, two more should be added which simply corre- 
spond to an interchange of the spaces X and Y (obtained by interchanging 
p and q). Although the interchangeability of these spaces was assumed in 
the derivation of inequalities (|76| ) and (^), clearly the sections A{x) will 
detect a difference. For example, there is no isometry which relates the sur- 
faces (10, 1,3) and (10,3, 1) (or their sections). On the other hand, the use 
of transformation (^) shows that the surface (10,3, 1) is isometric with the 
surface described by {n',p',q') = (5, 1,2) (where now p' < q'). However, not 
all interchanges of p and q produce a new Riemann surface. For example, 
{n,p,q) = (8,3,1) corresponds, by (|80D , to {n',p',q') = (8,1,3) (the values 
before the interchange). 

We summarize this discussion by our main Theorem: 

Theorem 3.6. Up to linear transformations, there are seven discrete Rie- 
mann surfaces [P\Q] generated by conformal maps of right triangles P and 
Q. Each of these surfaces is crystallographic; their properties, in particular 
the integers n, p and q which specify the angles of P and Q, are given in 
Table 1. 

Proof. The properties listed in Table 1 are simple consequences of general 
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results. Since 

G/A = (a,i?,r(7r,7r)), (90) 

we see that G/A is isomorphic to an abstract group generated by an element 
of order 2, a, and an element f, which has order 2\R\ = 2n, if n = \R\ is odd 
(and therefore r(7r, vr) ^ R) or n is even and just one of p and q is odd (since 
then the element of order two in R is either r(7r, 0) or r(0,'7r)). Otherwise 
(n even, p and q both odd), r(7r,7r) G R and f has order \R\ = n. These 
generators have the relation ara = and imply G/A ~ dn, the dihedral 
group of order 2n, or G/A ~ d2n- 

To compute the genus we use the orders of the vertex groups, (ffOD, (|8lD 



and (|82D, in formula (p^): 



2-2«=|G/A|(i- + i^-i). (91) 

The geometry of A is completely specified by the Gram matrix formed 
from its generators e^, /c = 1, . . . , 4: 

{Mki) = Ck ■ eu (92) 

where • is the standard inner product. If we consider each vector 
in and form the 4x4 matrix E whose rows are e^, then M = EE^^ 
and |A| = I det£^| = y/ det M. Using the LLL algorithm it was found that 
the generators of E and EE~^ (and hence A) could always be written as, 
respectively 

efc = 2(ae*2/ca^^gi2fc/3^) (gg) 
Cfc = 2 (a(e''" - l)e'''^", b{e'^^ - l)e'^^^) . (94) 

To help identify the lattice geometry it was sometimes necessary to define 
a new basis E' = SE, where 5* G SL (4, Z). The new Gram matrix is then 
given by M' = SMS^'^ . Details of this analysis, for the seven combinations 
of {n,p,q) in Table 1, are provided in the appendix. 

The only additional data needed in the density formula (|4^) is the triangle 
area \P\ = (1/4) sin27r(p/r;,). 
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Table 1: Properties of the seven discrete Riemann surfaces generated by 
conformal maps of right triangles 

The last column of Table 1 identifies which surfaces have doubly periodic 
sections. In general, since the kernel of the homomorphism vr^ : A — > 7rx(A) 
is given by the lattice Ay = A fl F, 

rk Ay = rk A - rk vrx (A) . (95) 

Surfaces with doubly periodic sections have rkAy = 2, while rkAy = 
corresponds to completely quasiperiodic sections. These are the only cases 
that occur, since (from (H) and (|9|)), 7rx(A) ~ ^[e^^^^P/")], the cyclotomic 
lattice with rank 0(n/ GCD {n,p)). □ 
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3.4 Piecewise flat surfaces and model sets 

Suppose a Riemann surface S is deformed into another surface, S, not nec- 
essarily representable locally by graphs of holomorphic functions. As long as 
the deformation preserves the isometry group and transversality with respect 
to Y, all the crystallographically relevant properties of the point set A{x) 
will be maintained in the corresponding deformed point set A{x). Kalugin's 
formula (|T7|), for example, makes this invariance explicit for the density. 

When a Riemann surface is generated by a triangular graph P\Q, defor- 
mations that preserve isometry and transversality are easily specified by the 
map defining the fundamental graph, f : P Q. We recall that in the 
Riemann surface, / is holomorphic and extends to a homeomorphism on the 
closure P. For the deformed surface S we continue to use the edge group 
G, defined by the geometry of the triangles P and Q, to form the orbit of 
the fundamental graph, but insist only that / is a homeomorphism. While 
still preserving isometry and transversality, we will even go a step further 
and modify the topology of S by defining S as the orbit under G of an open 
graph. Under these circumstances, when iS is a collection of disconnected 
components, we are free to relax the condition that / is a homeomorphism. 
In fact, we will primarily be interested in the case when / is a constant map, 
thereby making each piece of S fiat. 

Let P\y represent the graph of the constant map, f{P) = y. We will weigh 
the merits of various choices of y G F as giving optimal "approximations" 
of the map P\Q defined by the conformal map of triangles. By choosing 
y = Qi, a vertex of Q, we go the furthest in restoring partial connectedness 
to S. This is because the action of the vertex group Gi on P\Qi generates a 
flat polygon (possibly stellated) composed of triangles. Thus \Gi\ surface 
pieces will have been "aligned" by this choice of y. With the exception of the 
n = 5 and n = 10 surfaces, \Gi\ = {G-sl = 2n > IG2I, suggesting that either 
one of the singular vertices is a good choice for y. 

There is another criterion, however, that applies uniformly to all the 
surfaces and even distinguishes among the two singular vertices. A natural 
question to ask is: which value of y G Q "occurs with the highest frequency" 
in the graph P\Q7 To give this question a proper probabilistic interpretation, 
we suppose that x is sampled uniformly in P and ask for the probability that 
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\\f{x) — yo\\ < A, where / is the conformal map of triangles and A is the 
radius of a small disk about & Q. Since / is conformal, this condition (for 
A — > 0) is equivalent to ||a;— a;o|| < A/||/'(xo)||, where Xq = f~^{yo). Thus the 
probability of finding y in a neighborhood of yo is maximized by minimizing 
||/'(a;o)||. At a singular point ||/'(a;o)|| either vanishes or diverges, and in our 
case vanishes for xq Qi since we always have p < q. We therefore choose 
y = Qi for all of our surfaces. 

The result of flattening the seven Riemann surfaces in Table 1 by this 
prescription is particularly simple. In all cases V = Gi{P\Qi) is a regular 
{n/p)-gon covered p times (we restore connectedness to V by including edges 
incident to vertex 1). This means that if each point of A{x) is counted with 
multiplicity p, then A{x) and A{x) have the same density. Using Lemma 
we have 

S = A-V, A=(S), (96) 

if r(7r,7r) G R, otherwise, 

5 = A U At(2a,26)r(7r,7r) -P, A = (SS^^). (97) 

With the exception of the surface (5, 1, 2), r(7r, n) ■ V = V, and we have the 
simpler description: 

S = (S) V, (98) 

since 

(SS~^) U (SS-^)t(2a,26) = (S). (99) 

Figure 1 compares the point sets A{x) and A{x), given by the flattening 
process just described, of the Riemann surface (5, 1,2) of Table 1. Edges have 
been added to A{x) to aid in the visualization of three tile shapes: the boat, 
star, and jester's cap. There is a one-to-one correspondence between the two 
point sets, with most pairs having quite small separations. In any case, we 
are guaranteed the separation of corresponding points never exceeds 1, the 
diameter of triangle Q. A{x), the vertex set of a popular tiling model is 
a Delone set. A{x) fails to be a Delone set because triples of points appear 
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with arbitrarily short separations. That always three points coalesce in this 
way is a signature of the order of the branch points of ttx for this surface. 
A{x), on the other hand, has a "dynamical" advantage over A{x). Seen as 
atoms in a crystal or quasicrystal, the positions A{x) evolve continuously 
(in fact analytically) with x (viewed as a parameter), while the atoms in 
A{x) experience discontinuous "jumps", and for the most part never move 
at all. The singular loci of x G X for the two point sets have different 
dimensionalities: for A{x), 1 for A{x); the dynamics of A{x) is thus more 
regular also in this sense. To emphasize this point, we note that if 7(t) is 
almost any curve in X, then A{'~f{t)) is a regular homotopy (see Lemma ^.121) . 

The sets A{x) are examples of model sets, as defined by Moody 0, and 
therefore belong to the larger family of Meyer sets. That the set A{x) shown 
in Figure 1 can be organized into a finite set of tile shapes, for example, is a 
general property of model sets. Model sets obtained by flattening the other 
six Riemann surfaces of Table 1 are shown in Figure 3. In the three periodic 
cases, of course, even A{x) is a Meyer set and flattening is not necessary 
if that is our only goal. We flattened these surfaces because the sets A{x) 
are then particularly symmetric. In the three quasiperiodic cases the sets 
A{x) again organize themselves into tilings that have been discussed in the 
quasicrystal literature O, iMl O . 
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Figure 3: Point sets given by sections of flattened Riemann surfaces. Lines 
have been added to help identify tiles. 
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A Appendix 

For each of the entries in Table 1 we give below the corresponding Gram 
matrix for the generators (either (^) or (0)). From the transformed Gram 
matrices we see that the lattice geometries are in all cases simple root lat- 
tices {An, Dn, Zn) OT direct products. The lattices for (5,1,2) and (12,3,4) 
are obtained, respectively, from the lattices of (10,1,3) and (12,2,3) by an 
interchange of the spaces X and Y. 
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